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1 Introduction 



One of the major difficulties in developing quantum technologies, such as quantum 
computers JIJ, , are the different kinds of specifically quantum dynamical instabilities 
that can occur due to interactions between different degrees of freedom and resonant 
interaction with the external fields. These instabilities differ from dynamical insta- 
bilities in classical systems, which are usually connected with strong dependence of 
trajectories on the initial conditions and on the values of parameters. Small variations 
of initial conditions or parameters lead to large deviations in time of the corresponding 
trajectories. If the speed of this deviation is exponential, the system becomes chaotic, 
and the appropriate methods of description are statistical rather than deterministic. 
But for quantum systems, the notion of a trajectory is not well defined. This is one 
of the main reasons why most of the well-developed methods for stability analysis can 
not be directly applied to quantum systems. Moreover, as was first shown theoretically 
by Berman and Zaslavsky || |J (see also ||), even in a "deep" quasiclassical region, 
classically chaotic systems can have quantum dynamics that is very different from the 
corresponding classical dynamics. 

Another important phenomenon which takes place in quantum systems which are 
classically chaotic is quantum nonlinear resonance (QNR), which was first introduced 
and investigated theoretically by Berman and Zaslavsky || . QNRs are quantum man- 
ifestation of nonlinear resonances which play very important role in classically chaotic 
systems 0-103 • Interactions of QNRs are known to be intimately connected to the 
transition to quantum chaos ||11||-||17||. In the simplest situations, QNRs occur when a 
bound quantum system whose energy levels are not equally spaced is driven by a res- 
onant perturbation. A QNR is characterized by two main parameters: the number of 
quasi-energy levels, 5n, which are "trapped" in the potential well of the resonance, and 
the characteristic frequency of slow phase oscillations, VL ph . Isolated QNRs imply stable 
quantum dynamics; overlapping QNRs cause a transition to quantum chaos. QNRs are 
very general phenomena in non-integrable quantum systems, and can be thought of as 
"quasi-particles" of quantum chaos (for more details, see chapter 9 of Reichl's recent 
book |l8j, which is devoted to the transition to quantum chaos caused by interaction of 
QNRs). Until now, QNR effects have been experimentally investigated using Rydberg 
atoms in a resonant micro- wave field |T9[]. Understanding the instabilities connected 
with overlapping QNRs is important for fundamental problems related to the transi- 
tion to quantum chaos, and for the design of experimental devices (such as quantum 
computers based on ion traps) in which these instabilities may cause significant prob- 
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lems. To study the characteristic parameters of both isolated QNR's and the problems 
related to interaction of QNRs it is important to choose a model which (a) involves 
regulated (and relatively small) number of interacting QNRs; (b) can be implemented 
experimentally in quantum and quasiclassical regions of parameters. 

In this paper we introduce a quantum model which is convenient for investigation of 
quantum dynamical instabilities and transition to quantum chaos based on overlapping 
of QNRs. The model consists of a single ion confined in a radio-frequency Paul trap 
and which interacts with a resonant laser field. In the classical limit, this model reduces 
to the well-known model of a linear oscillator interacting with a plane electromagnetic 
wave, and was investigated in |20) (see also references therein). The main advantage of 
our model is that the number of interacting QNRs can be regulated, for example, by 
varying the intensity of the laser radiation, which is difficult to achieve in other models 
based on the kicked interaction ||24||-||28||. 

Devices based on trapped ions have been used to investigate experimentally fun- 



damental aspects of quantum mechanics [29, 30], as well as for important technological 
applications such as optical frequency standards |yj and quantum computing |], |32| . 
Ions are confined by a combination of a rotating quadrupole potential (induced by 
the rod electrodes) and a weak electrostatic potential (induced by the conical endcap 
electrodes). The ions, once trapped, can be cooled by standard Doppler cooling and 
by a optical pumping method ("sideband cooling"), which can cool multiple ions down 
to the quantum mechanical ground state of the trapping potential. In an ion trap 
quantum computer, information can be stored in the internal quantum states of the 
ions (which constitute the quantum bits, or "qubits" of the computer), and, using ul- 
tra narrow bandwidth lasers, quantum gate operations can be realized between pairs 
of qubits using quantum states of the collective motion of the ions in the harmonic 
confining potential as a quantum information bus [33]. As such devices are specifically 
designed to investigate experimentally the preparation, evolution and measurement of 
quantum systems with large dimension Hilbert spaces, the linear ion trap is an ideal 
apparatus to investigate the problems of quantum dynamical stability, the transition to 
quantum chaos, and the spectroscopy of quantum nonlinear resonances. In this paper, 
we present the main elements of the derivation of our model - a quantum linear oscil- 
lator driven by a monochromatic wave, and the preliminary analytical and numerical 
results on the classical and quantum dynamics in different regions of parameters. 

The paper is organized as follows. In section 2, the theory of how a trapped ion 
can be driven by laser fields in the manner of the harmonically driven oscillator is 
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described in detail. It should be noted in particular how similar the arrangement and 
laser requirements are to those employed in ion trap quantum computer experiments. 
In section 3, the classical theory of the harmonically driven oscillator is discussed; the 
quantum theory is described in section 4. The connection of this system with the solid- 
state Anderson localization model is described in section 5. The results of numerical 
simulations are presented in section 6. We conclude this paper with a brief discussion 
of the possibilities for experimental verifications. 



2 Raman interactions of lasers and trapped ions 

A single ion confined in a linear radio-frequency (rf) trap may be described by an 
effective Hamiltonian given by the formula 



H = —p 2 + -mu 2 x 2 + H I} (2.1) 
2m 2 

where m is the mass of the ion, x (p) is the position (momentum) operator for the 
ion and uo is the angular trapping frequency. We are only considering motion of the 
ion along one direction, namely the axis of weak confinement of the trap; the ion is 
strongly confined along the other two directions, transverse to the axis and so we will 
assume that the motion in those directions can be neglected (see Fig.l). 

We will employ the interference of two laser beams acting on the ion to realize 
experimentally our desired interaction Hamiltonian Hj . Such Raman interactions be- 
tween lasers and ions are a standard technique, and are described in detail elsewhere 
34j|. The ion, confined in the harmonic trapping potential, will have many quantum 



levels associated with both the internal (atomic) variables and the external (motional) 
degrees of freedom. We will confine our attention two manifolds of such states, sepa- 
rated in energy by an appreciable amount (see Fig.2). What we have in mind is a lower 
manifold consisting of the magnetic sublevels of our ion, each level having a series of 
sidebands associated with excitation of quanta of the external harmonic oscillations; 
the upper manifold would then be the sublevels of an excited state of the ion, with 
similar sidebands. The lasers with which the ion is interacting will be detuned from 
the optical transition between the upper and lower levels, so that there is a negligible 
probability of any of these levels becoming excited: the lasers only cause redistribution 
of population amongst the lower manifold of levels. The upper levels may then be 
adiabatically eliminated from the problem, and one can therefore show that the matrix 
elements of the effective interaction Hamiltonian for the lower manifold is given by the 
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formula [|3l 



(M\Hj\N) 



EE 




exp [i {w a -Wat)t], 



(2.2) 



L a,a' 



where the sum involving u is over all of the upper manifold levels and the two sums 
involving a and a' are over all of the applied laser fields, the Rabi frequency of the 
a-th laser being defined by 



In Eq. ( p.3|) , di is the z-th component of the dipole moment operator ( i = (1, 2, 3), 
standing for the three Cartesian components of a vector, and summation over repeated 
indices being implied), e\°^ is the z-th component of the electric field from the a-th 
laser (which is a function of the ion's position operator, r), Tilom is the energy of the 
M-th lower manifold level, Hu>l is the energy of the L-th upper manifold level, and u a 
is the angular frequency of the a-th laser. 

To proceed, we will make the distinction between internal and external degrees 
of freedom. We can form a basis set for the Hilbert space from a tensor product of a 
set of internal quantum levels with a basis set for the external degrees of freedom (for 
example the Fock states of the Harmonic oscillator). The set of internal states will 
be divided between the upper and lower manifolds. Thus we will make the following 
substitution: 



where \m), \n) and |/) are members of the basis states for the motion degrees of freedom, 
| A) is a member of the upper internal manifold and \fj) and \v) are members of the 
lower internal manifold. In this notation, the matrix elements of the Hamiltonian Eq. 
Q become: 



(2.3) 



\M) 
\N) 



|A)|0 

\n)\m) 

\v)\n), 



(2.4) 
(2.5) 
(2.6) 



{n\{m\H!\n)\v) 




^\di\\){\\dj\v) 



ATi(uj\ -lo v — uj a i +Ut — u n ) 



(m\Ej a \v)\l)(l\Ef>(r)\n) exp [i (u a - u a ,) t] . 



(2.7) 



The average detuning of the A-th upper manifold level is be defined to be 



= LO\ — lo v — u, 



(2.8) 
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where %uj v is the average energy of the lower manifold and ui a is the average of the 
laser frequencies. We will assume that, in the denominator of the fraction appearing 
in Eq. ( |2.7|) , we can make the following approximation: 

u\ - u u - to a > +uji-uj n K, A X . (2.9) 

If we use the completeness property of the external basis states (i.e. J2i 10 (^1 — Ii where 
I is the identity operator), then we obtain the following formula for the Hamiltonian 
operator for the lower manifold states: 

= MR (2-10) 

where 

(f , 4) = "E W^W&M e, (r, t) E- (r, t) , (2.11) 
and the total laser field Ei is the sum of the different laser components: 

Ei (r, t) = J2 E i a \r) exp (icu a t) . (2.12) 

a 

2.1 Two level systems 

Let us now assume that there are only two internal levels in the lower manifold, which 
we will denote |1) and |2). As will be discussed below, this is a reasonable assumption to 
make for the atomic systems we have in mind. Also, we introduce a special coordinate 
system: the two internal levels are split by a magnetic field acting along the Z axis, 
which is the axis of quantization for the internal levels of our ion. The other two axes 
are the X and Y axes. These axes do not necessarily coincide with the x, y, z directions 
introduced to describe the motion of the ion in the trap. In this case it is convenient 
to use the Pauli operators for the system: 

a x = |1)(2| + |2)(1|, (2.13) 
a 2 = i(|l)(2|-|2)(l|), (2.14) 
a 3 = |2)(2|-|1)(1|. (2.15) 

Using these operators the Hamiltonian can be written as follows: 

H I = ho(r,t)I + h i (r,t)a i , (2.16) 
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where 7 is the identity operator (|1)(1| + |2)(2|) and 



h (r,t) = (r,t) + «2, 2 (r,t)], (2.17) 

hi(r,t) = -[K 1>2 (r,t)+K 2>1 (r,t)}, (2.18) 

h 2 (i,t) = -[K 1)2 (r,t)-K 2il (r,t)], (2.19) 

fc 3 (r,t) = ~[K 2 , 2 (r,*)-«i,i (*,*)]■ (2.20) 



For the special case that the lower manifold of internal states consists of two 
magnetic sublevels of the 2 Si/2 ground state of an alkali-like ion, and the upper manifold 
is the two sublevels of the 2 P\/2 excited state whose Zeeman splitting is not too big, 
the atomic matrix elements appearing in Eq. (|2.11| ) can be calculated in closed form. 
As a result the components of h appearing in Eq. ( |2.16| ) reduce to the following simple 
form (see appendix A): 

h (r,t) = - X |E(r,t)| 2 

/H(f,t) = 2 X lra{E z (i,t)E* Y (v,t)} 

h*(r,t) = 2 X lm{E x (r,t)E* z (r,t)} 

h 3 (r,t) = 2 X lm{E x (r,t)E* Y (r,t)}, (2.21) 

where Im{...} is the imaginary part of the quantity in curly brackets and x = 
Ane /4:kQA (k and A being, respectively, the wavenumber and the Einstein A coeffi- 
cient for the transition between the upper and lower manifolds, A the laser detuning 
and e the permittivity of free space). 

The quantity proportional to ho (r, t) in Eq. fl2.16| ) represents a dynamical effect 
of the laser fields on the ion which does not cause any effect on its internal degrees 
of freedom; the term proportional to hs (r,t) represents a A.C. Stark shift of the two 
internals levels; the terms proportional to hi (r, t) and /i 2 (r, t) represent transitions 
between the two levels of the lower manifold. If we make the requirement that the 
lasers are plane-polarized along the axis of quantization Z, then it is clear from the 
above formulas that h\ = hi — h% = and only the first term involving ho has any 
effect. 

Let us assume that two laser beams, designated the pump (p) and the Stokes (s) 
beams are present, both plane polarized in the Z-direction, i.e. 

E x {r,t) = 0, 
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E Y (r,t) = 0, 

E z (r,t) = E ip) exp [-i (k p • r - u p t)) + E (s) exp [-i (k s ■ f - cu s t)] . (2.22) 
The interaction Hamiltonian in this case is given by 
Hj = x {\e^\ 2 + \E^\ 2 + 2 \E^E^*\ cos [(k p - k s ) • r - (u p - u a ) t + 0]} , (2.23) 

where <p — Arg E^*^ is the phase difference between the two lasers. The constant 



terms involving 



E(p) 



and 



E (s) 



have no effect on the evolution, and so will be 
neglected. Thus the full Hamiltonian, including the effect of the harmonic evolution of 
the ion along the weak axis of the trap (but excluding the internal free evolution) is 



~2 2^2 

p muj x e „ , 

H = — H h — cos(A;x - fit). 

2m 2 k J 



(2.24) 



where Q = u p — uo s . The parameters e and k, which will feature prominently in what 
follows, are given by 



e 
k 



AlT€ok 

~2kJK 
(kp — k s 



_g(p)_g( s > 



(2.25) 
(2.26) 



where e x is the unit vector along the x-axis, i.e. the axis of weak confinement in the 
trap. 

3 Classical Limit 



In the classical limit (p — > p, x — > x) , the Hamiltonian ( p.24|) takes the form, 

2 2 2 

p muj x e ,, _ . 

ti = 1 h — cosffcx - \lt), 

2m 2 k 

The classical equations of motion in (p, x) variables are 



V 



dH 



—muj 2 x + esin(kx — Qt), x 



dH p 



(3.1) 



Equations 



dx dp m 

lead to the following second order non-linear differential equation 



(3.2) 



X + U) X 



m 



sin (kx — fit). 



(3.3) 



S 



3.1 Dynamics near resonances 

Assume that the driving frequency is close to a resonance, i.e. 



Nu « fi, (3.4) 

where N is an integer. In this case, it is convenient to describe a classical dynamics 
using the "action-angle" variables f2(| [21], ^], (J, <p), which are related to the variables 
(p,x) by the canonical transformation (see Appendix B), 



2Ni fip + nt\ , — - — /<p + nt\ 

x = \ cos — — — , p = — v 2N Iiumsml — — — . (3.5) 

V urn \ N ) \ N 1 

In the variables {I,<p), the Hamiltonian ( |3.1|) takes the form, 



H = (Nu + r cos ( k\^— cos $ - fit V (3.6) 

k \ V urn J 

where 

* = (3.7) 



The second term in (|3.6| ) can be expanded as a series of Bessel functions J n (z) 



/ 2NI \ 00 / 2NI\ 
cos W cos$-fit = ^ Jn W cos[n($ + vr/2) -fit]. (3.8) 



We have for the phase, n($ + vr/2) — fit, 

n($ + 7r/2) - fit = n-^: + — + — fit. (3.9) 

Thus, the classical Hamiltonian ( ^.6[ ) can be represented in the form where the unper- 
turbed part and the perturbation are explicitly separated, 

H = H + H int , (3.10) 

where 



H = (Nu-n)I+^J N (k\^^)cos(cp + Nir/2), (3.11) 

k \ V rnuj / 



. ^ , 2NI\ n nn n-N„\ 

H*nt = rE^lW— cos(-^ + T + ^ r fit . (3.12) 
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In the (J, if) variables, the classical equations of motion are, 
dH 



dip 



dH 
~dT 



k 



J N (z) sm((f + nN/2) + 



e 



n 



— J2nJ n (z)sm[-ip + — + 



nn n — N 



kN . 

NiU -Q + S\ 



N 



N 



2mul 



J' N (z) cos((p + irN/2) + 



N 



2mujl 



rm n — 



N' 



N 



where, 



z = k\ 



I2NI 



(3.13) 
(3.14) 

(3.15) 
(3.16) 

(3.17) 



is a dimensionless variable. Equations ( ft.14 ) and (|3.16| ) are convenient when analyzing 
the classical dynamics in the vicinity of the resonance ( |3.4j ). This case corresponds to 
small values of e, 

e < e cr . (3.18) 

Usually, the critical parameter e c in (3.18) should be found numerically (see Sec. ||). 
Under the condition ( |3.18| ), classical dynamics can be approximately described by the 
Hamiltonian H ( |3.11|) . The corresponding approximate equations of motion follow 
from Eqs (|3~T4T ) and ( PlBj ), 

• dH e 



0, 



k 



J N (z) sm(cp + irN/2), 



— = Nu J -Q + e ] 



N 



2mul 



J' N {z) cos(v9 + TlN/2). 



(3.19) 
(3.20) 



In the general case (for e large), it is more convenient to use the following exact 
equations written in (J, ip) variables, 



e 2NI 



N V moo 
Nuj-Q-e\ 



sin 



Am 



I2NI fp + Qt 
cos 



N 



Qt 



sin 



V + Qt 



N 



N 



2mul 



sin 



I2NI fip + Qt 
cos 



N 



Qt 



cos 



'tp + Qt^ 
N 



(3.21) 



(3.22) 
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3.2 Dimensionless variables 



To describe both the classical and quantum dynamics, it is convenient to introduce the 
following dimensionless variables, 



n 



T = U)t, 

H 



o 



n 



Hint 



3 ■ -int j_ 

nuj noo 
and the dimensionless parameters, 



n 



n 



i 

" n' 

H 



(3.23) 



V 



fi = — , 8 = N — fjt. 



(3.24) 



mu^ 2muj ' " u 

The parameter r] is the Lamb-Dicke parameter used in the theory of ion traps to 
quantify the strength of confinement. It is related to effective Planck constant by the 
formula % = 2r) 2 . 

3.3 Isolated Nonlinear Resonance 



Using (|3.23| ) and ( |3.24|) , we have from ( |3.19| ) and (|3.20|) the approximate dimensionless 
equations of motion in the vicinity of the resonance (|3.4|) , 

d£ dHo e _ , . . . _ T/r> . 

Tr = - a7=v* Wsm( ^ + " JV/2) ' 

dip dH 
~dr~ ~ ~W 



e N 

^ + Y\I- J ^ Z ) cos(^ + 7riV/2), 



(3.25) 
(3.26) 



(3.27) 



where z = 2r\\/Ni and the dimensionless resonant Hamiltonian is, 

Ho = £8 + t^Jn(z) cos(^ + ttJV/2). 

The classical dynamics corresponding the Hamiltonian ( 3.27]) we shall call "Nonlinear 
Resonance" (NR). 

To estimate the region of parameters of validity of Eqs (|3.25 ) and ( |3.26|) , their 
solutions should be compared with the solutions of exact equations ( |3.21| ) and ( |3.22| ). 
Equations ( |3.21|) and fl3.22|) in the dimensionless variables have the form, 



dr 



rj V N 
e 

2q 



sm 



z cos 



dp) , e IN . 
—r- = 8 — — \ I — sm 
dr 



z cos 



s N )' 

'ip + [IT 



fir 



sm 



N 



jJT 



, at y 

' \p + fir' 



COS 



N 



(3.28) 
(3.29) 
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Equations ( |3.28| ) and Q3.29D are derived from the exact dimensionless Hamiltonian, 



H = £5 + — cos 
2r] z 



z cos 



N 



In the dimensionless variables ( |3.23| ) and ( |3-24 ) , Eq. 



jJLT . 

takes the form, 



(3.30) 



d 2 ^ 

— + £ = esin(£-/ir). 



(3.31) 



4 Quantum Equations of Motion 



In the dimensionless notation (|3.23j ) and ( |3.24|) , the quantum Hamiltonian ( |2.24|) takes 
the following form in the coordinate representation, 



— = n = — 

hu 2rf 



d 



- 2rj —-r + — + e cos(£ - fir) 



ae ' 2 

The Schrodinger equation for the Hamiltonian (^4.1| ) is, 

2^(C,r) 



(4.1) 



i2rf 



dr 



Hlo + e cos(£ - jir) 



(4.2) 



where fiio is the Hamiltonian of a linear oscillator, 

d 2 £ 2 

For ho we have the well-known eigenvalue problem, 

H LO \n) = 2 V 2 (n + l/2)\n), 



(4.3) 



where, 



\n) = 0n(O 



1/2 



H n (Z/V2 V )e 



2 /V 



(4.4) 



(4.5) 



j]2 n n\V2n_ 

where H n (y) is a Hermite polynomial. Although these eigenfunction may appear some- 
what unfamiliar because of the use of dimensionless variables, they are in fact the 
standard eigenfunctions of an unperturbed harmonic oscillator (i.e. the Hamiltonian 
given by Eq. (2.1) with Hj = 0). The normalization condition for the eigenfunction 

0n(O IS, 

*n(£)0m(£R = S n,m- (4.6) 
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To describe the quantum dynamics we represent the wave function in (|4. 2|) in the 



form, 



^,r) = ]>>n(r)0 n (O- 

n=0 



(4.7) 



^From (|4.2|) we have the equations for the complex amplitudes, c n (r) 



(m + l/2)c m (r) + — J2( m \ cos (£ ~ ^r)\n)c n {r) 

I n=0 



■ dCrJj) 

dr 



= (m + l/2)c m (r) + — £ ( e-^F^fo) + e^F^J^J c n (r). (4.8) 

In Eqs (f4.8|), F m>n (r/) is the matrix element, 

F mw fa) = (m|e*|n) = -p= / H m (u)H n (u)e- u +l2vu du. (4.9) 

Equations ( |4.8| ) are used below, in Sec. |5], for numerical calculation of the quantum 
dynamics of the system. 



5 Results of Numerical Calculations 

In this section, we present results of numerical simulations of classical and quantum 
dynamics of systems whose Hamiltonians have the form given by Eq. ( p. 24 ). We 
have used a set of parameters which will allow easy experimental verification of our 
predictions using the type of ion trap apparatus currently being used to investigate 
quantum computation. 

If we use the geometry for the pump and Stokes lasers shown in Fig.l, the pa- 



rameter k defined by Eq. ( |2.26|) is given by k = cos6(k p + k s ) pa 2/cqcos^. The laser 



field strengths 



and 



can be related to the power in the pump and the Stokes 
beams respectively. It is usual to generate one of these beams (the Stokes, say) by 
frequency modulation of the pump beam, so that the beam parameters will be similar 



for them both. The power in the pump beam is given by the formula (Ref. [P5], p. 
488) 



P 



ce ?r 



-Wn 



E (p) 



(5.10) 



where w is the laser spot size (i.e. the 1/e 2 radius of the intensity distribution). If 



we substitute this result into Eq. ( 2.25 ), we obtain the following expression for the 
dimensionless parameter e: 

ek 
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where s = / is a dimensionless parameter of order unity which can be 

controlled experimentally. 

Singly ionized calcium (m = 6.64 x 10~ 26 kg) is common ion in use by several 
groups worldwide for ion-trap quantum research (see Fig. 3, for energy levels of this 
ion). For the 2 S 1 / 2 — 2 Pi/2 transition in this ion (wavelength Ao = 397nm) the Einstein A 
coefficient is |J3| A = 1.30 x 10 8 s _1 and the wavenumber is ko = 2ti/\q = 1.58x10 m -1 . 



If we assume a laser power of lOmW, a spot size wq = 20/zm a trapping frequency 
u = 2tt x 500kHz, and a detuning A = 2vr x 1.0GHz, then e = 1333s cos 2 9. Thus 
by varying the experimental free parameters s and 9 one can achieve a large dynamic 
range for the dimensionless driving force e. The Lamb-Dicke parameter rj for these 
parameters is 0.502 cos 9. 

The driving frequency Q is the difference between the pump and the Stokes fre- 
quencies. As mentioned above, these two beams will be realized by splitting one parent 
beam using a beam splitter and then frequency modulating one of the resultant beam 
using either an acousto-optic or an electro-optic modulator. In this manner splittings 
as high as Q = 2tt x 100MHz can be achieved without too much difficulty, so that the 
dimensionless parameter N can be in the range to 200. 

5.1 Simulation of classical dynamics 

In Figs 4a and 4b, the classical phase space (Poincare section) is shown in (£ , <i£ / dr) 
plane, for seven initial conditions, and for different values of the dimensionless driving 
force, e, which characterizes the intensity of the laser beams, and is defined in ( |5 . 1 1| ) . 
To derive these results, the Eq. ( |3.31|) was solved numerically (for N = A and 5 = 10~ 2 ). 
One can see from Fig. 4a, that for small values of e (e < 2) the classical dynamics is 
regular in some regions of the phase space. Note, that even at these values of e there 
exist relatively large regions in the phase space with the chaotic component. When 
the interaction parameter e increases, the regions with the regular classical dynamics 
become smaller. As one can see from Fig. 4a, even at e = 4, the dynamics in the region 
of the phase space corresponding to the "classical ground state" (CGS) (the vicinity of 
the point (0, 0)) remains regular. At larger values of e (e ~ 8) the CGS becomes chaotic. 
This transition is demonstrated in Fig. 4b. Figs 5 and 6, show the time evolution of 
the dimensionless classical coordinate of the ion, £(t). The maximum dimensionless 
time of simulations, r max = 100, corresponds to the real time-scale t max = 31.8/is 



14 



(u = 2tt x 500kHz). The initial conditions for the dynamics shown in Figs 5 and 6 
correspond to the CGS. As one can see from Fig. 5, the chaotic component appears at 
e ?s 8, and is well resolved for e > 10 (see Fig. 6). The frequency Fourier transform for 
the dynamics shown in Figs 5 and 6, is presented in Figs 7 and 8. It is well-known, that 
the transition to the dynamical chaos in classical dynamical systems is accompanied by 
the transition in the frequency Fourier spectrum. The regular dynamics corresponds to 
the discrete frequency spectrum, the chaotic dynamics corresponds to the continuous 
frequency spectrum. This characteristic modification of the frequency spectrum is 
demonstrated in Figs 7 and 8. For small values of the dimensionless driving force e, 
one can see only some quasi-discrete lines in the Fourier spectrum. When e increases, 
the frequency spectrum transforms to the continuous one. 

5.2 Simulation of quantum dynamics 

To simulate a quantum dynamics the following parameters were chosen: N — 4, 
5 = 10~ 2 , rj = 0.45. For the initial conditions we used the ground state of the un- 
perturbed quantum linear oscillator: c (0) = 1, c n (0) = for n > 0. Fig. 9, shows the 
time evolution of the quantum probabilities P n (r) = |c n (r)| 2 (n=0, 1, 2, 3, 4 and 5) 
for relatively small values of e =0, 0.5, 1, 1.5 and 2. These values of e correspond to 
the regular classical dynamics which starts from the CGS. Time evolution of quantum 
probabilities P n (r) (n — 0, 1, 2 and 3) for bigger values of e is shown in Fig. 10. 
Because the value of the parameter e ps 8 corresponds to the classical chaotic dynam- 
ics for the initially populated CGS, the curves (4) in Fig. 10 describe the quantum 
chaotic motion. Fig. 11, shows the dynamics of the probability function Po(t) for the 
larger time interval: r G [0,30]. In the real time this correspond to: t G [0,9.54]/xs. 
As one can see from Fig. 11, for e > 7.5 the dynamics of the probability function 
Po{t) becomes rather complicated, and correspond to the classical chaotic motion. We 
hope that this complicated dynamics of the probability functions P n { T ) can be mea- 
sured directly in the experiments with trapped ion. Fig. 12, represents the results of 
the numerical simulations of the frequency Fourier spectrum, Po(v), of the quantum 
probability function Po{t). As one can see from Fig. 12, the characteristic qualitative 
modification of the frequency spectrum P^v) starts from e > 7.5. This modification 
of the frequency Fourier spectrum can also be measured experimentally. In Fig. 13, 
we show the results of numerical simulation of the dynamical evolution of the average 
value: (£ 2 ) = ( l I / (^, r)|£ 2 |\I/(£, r)), where the wave function \I/(£,t) is defined in (4.7). 
This dynamical characteristic is very important for understanding the conditions of 
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stability of the system under consideration, as it describes the amplitude of the ion's 
oscillations in a trap due to the influence of the resonant laser fields. Experimental 
measuring the time dependence of this amplitude is important for characterizing the 
regular and chaotic dynamics of an ion. The frequency spectrum of the amplitude (£ 2 ) 
is shown in Fig. 14 for different values of the perturbation parameter e. As one can see 
from Fig. 14, at e > 7.5, the frequency spectrum qualitatively modifies, and includes 
many harmonics. In Fig. 15, we show the time evolution of two functions: the unper- 
turbed Hamiltonian of quantum linear oscillator (4.3), (TIlo), and the amplitude (£ 2 ), 
for the values of the perturbation parameter e = 0; 05; 2. The corresponding classical 
dynamics is regular in this case. The time evolution of the same functions is shown 
in Fig. 16 for larger values of the parameter e = 3; 5; 7.5; 8. The curves (3) and (4) 
correspond to the classically chaotic regime of motion. The maximum simulation time 
at Fig. 16 is: t max = 4.77/^s. 

6 Conclusion 

In this paper, we introduced a quantum model which describes a transition from a regu- 
lar dynamics to quantum chaos, for a single ion in a linear ion trap. The configuration 
of the resonant laser fields allowed us to represent the model in a "standard" form. 
Namely, our model formally describes a quantum linear oscillator interacting with a 
one-dimensional plane wave (see the Hamiltonian (2.29)). The classical version of this 
model is very important and useful in the theory of the dynamical chaos 0. This 
model differs significantly from the "usual" nonlinear models (as, for example, a "stan- 
dard map" |8], [IIJ) considered in the theoretical works on dynamical chaos. Namely, the 
model described by the Hamiltonian (2.19) includes "nonlinearity" and "perturbation" 
in the same term. In this sense, this model is degenerated, and possesses many unusual 
properties ||. So, quantum analysis of this model both theoretically and experimen- 
tally will be of significant importance for understanding complicated dynamics in this 
system, and for applications in different devices based on the trapped ions, including 
quantum computer. The numerical results presented in this paper describe both regu- 
lar and chaotic quantum dynamics which starts from the initial population of the linear 
oscillator's ground state. We believe that this regime of motion should be investigated 
in the experiments with the trapped single ion. These experiments will allow one to 
establish qualitative and quantitative correspondence between the quantum dynamics 
described by the Hamiltonian (2.29) and the real system "a single trapped ion + res- 
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onant laser fields". Further theoretical analysis and experiments should include both 
pure quantum and quasiclassical regions of initial population, and the parameters de- 
scribing regular and chaotic regimes in both these regions. These investigations are 
now in progress. 
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Appendix A: Derivation of Eq. (|272T]) 



The dipole matrix elements appearing in Eq. (|2.11 ) can be written in terms of the 



Einstein A coefficient between the upper and lower manifolds as follows: 



Mim = f A ^ +1) t J " 1 Me'* (A.l) 

where c is the speed of light, a is the fine structure constant, u is the angular frequency 
of the transition between the upper and lower manifolds, (J\,m\) are the magnetic 
quantum numbers for the upper manifold state |A) ( J^, m^) are the magnetic quantum 
numbers for the lower manifold state the term contain six quantities in brackets is 
the Wigner 3-j symbol, and ef is the i-th component of the g-th normalized spherical 
basis vector, viz. 

e« = -i=(l,-i,0), (A.2) 

e<°> = (0,0,1), (A.3) 
eC- 1 ) = ^(M,0). (A.4) 

(A.5) 



Substituting Eq. ( A.l ) into Eq. ( [2.11 ), and assuming that the Zeeman splitting in the 
upper manifold is small compared to the overall detuning, so that A a ~ A (independent 
of A) , then we find that the quantities are given by the following formula: 

*w (r, t) = ^^Ay (A*, ») Ei (f , t) E* (r, t) , (A.6) 
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where the tensor A# (fi, v) is given by 



A,,- fav) = (2 J x + 1) £ E 

mx=-J\ q,q'=-l 



1 J. 



A 

-m„ q m x 



Ju 



1 J. 



A ] e (<l) e («')* 



(A.7) 

If we assume that the lower manifold is the 2 S1/2 ground state of an alkali ion, the 
two states being denoted |1) = 2 Si/2,rn = —1/2 and |2) = 2 Si/2,rn = 1/2, and that 
the upper manifold is the 2 P\/2 state, then these tensors can be found in closed form: 



MM) 



A?-(2,2) = i| ] 

' • 1 



1 -i 
1 



A.- (1,2) = -A*. (2,1) 




(A.8) 



Therefore the cross products involving the electric field components with these tensors 
can be written as follows: 





(1,1) EiE* = 


1 r 

3 


|E| 2 + 2Im{E x ^}" 




(A.9) 


A* 


(2,2) E t E* = 


1 r 

3 


|E| 2 -2Im{E x E$r} 


> 


(A.10) 




(1,2) E,E* = 


i [E z (E* x + iE* Y ) - E* z (E x + iE Y )\ , 


(A.ll) 


Aii 


(2,1) EiE* = 


h 


-E z (E* x - LE^) + E* z (E x - iE Y )\ . 


(A.12) 



If these results are substituted into the definition of h and hi, Eqs (|2.17|) - (|2.20|) , 
one obtains Eq. ( |2.21| ). Similar results, with slightly different numerical factors, are 
obtained if the upper manifold is a 2 P3/2 state. 



Appendix B: Canonical Transform to "action-angle" 
variables 



The theory of canonical transforms in classical mechanics is described in detail in 
Ref. [0, §45. We want to transform from a set of variable i,p to a new set (p, I, 
where (p plays the role of position coordinate and / the role of momentum. Such 
transforms are specified by a generating function, F(x,ip,t). Then variables p and I 
and the Hamiltonian in the new coordinate system are are related to F by the following 
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formulas: 

dF 

p - ^ (B.1) 

dF 

H( V ,I) = H(x,p) + ^ (B.3) 
For the canonical transform used in section 3.1, the generating function is given 

by 

t) = -—x 2 tan (^^J • (B.4) 



Substituting, we find: 



p = — mcuxtanl — — — J (B.5) 



2iV V AT 

= H(x,p)-QL (B.8) 
Equation (|B.5|) implies that 



z = W cos r AT . (B.9) 



muj V N 

On substitution of this last formula into Eq. (|B.6|) we obtain 



P 



-V2Nmuj I sin C^~jf^ ) • ( B - 10 ) 
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Figure captions 



Fig. 1: A schematic diagram of an ion in a linear trap to illustrate the notation and 
configurations of the laser fields. 

Fig. 2: A schematic illustration of the energy levels of a trapped ion. 

Fig. 3: Energy levels of Ca + ion. Wavelengths and radiative lifetimes are shown. See 
for references. 

Fig. 4: Classical phase space (Poincare section). Trajectories with different initial 
conditions are shown. The values of e are indicated in the figure: (a) e = 2,2.5,3,4; 
(b) e = 5, 8, 10, 20; r] = 0.45; N = A, S = 10~ 2 . 

Fig. 5: Time evolution of the dimensionless classical amplitude of oscillations £(r) 
for different values of e = 0, 0.5, 1, 2, 5, 8. The maximum time of simulation it t max = 
31.8/is. Initial condition: (£ , d£ / dr) = (0,0). 

Fig. 6: Chaotic dynamics of the dimensionless classical amplitude of oscillations £(r) 
for the values of e = 10, 20, 30, 40. The maximum time of simulation it t max = 31.8/xs. 
Initial condition: (£ , d£ / dr) = (0,0). 

Fig. 7: Frequency Fourier spectrum of the classical amplitude £(r) for e=0,0.5,l,2,5 
and 8. Transition to chaos appears at e « 8. Initial condition: (£,d£/dr) = (0,0). 

Fig. 8: Frequency Fourier spectrum of the chaotic dynamics of the classical amplitude 
£(r) for e =10,20,30 and 40. Transition to chaos appears at e ~8. Initial condition: 
(£,de/dr) = (0,0). 

Fig. 9: Dynamics of the quantum probabilities: -R«(t) = |c n (r)| 2 , for n =0,1,2,3,4 and 
5, and for e = 0, 0.5, 1, 1.5 and 2. Initial condition: co(0) = 1, c n (0) = for (n > 0). 

Fig. 10: Dynamics of the quantum probabilities: P n {T~) = l c n(T)| 2 , for n =0, 1, 2 and 
3, and for e =3, 5 and 7.5. Initial condition: co(0) = 1, c n (0) = for (n > 0). 

Fig. 11: Dynamics of the quantum probability: Po(r) = \c (t)\ 2 , for e =1, 5, 7.5 and 
8. Initial condition: c (0) = 1, c n (0) = for {n > 0). Transition to quantum chaos 
corresponds to e ~ 8. Initial conditions as in Fig. 10. 

Fig. 12: Frequency Fourier spectrum of the quantum probability: Pq{j) = \c (r)\ 2 
for e =1, 5, 7.5 and 8. Transition to chaos appears at e ~ 8. Initial conditions as in 
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Fig. 10. 

Fig. 13: Dynamical evolution of the average value: (£ 2 ) = r)|£ 2 |\I/(£, r)), where 
the wave function r) is defined in (4.7). e =3, 5, 7.5 and 8. Transition to chaos 
appears at e ps 8. Initial conditions as in Fig. 10. 

Fig. 14: Frequency Fourier spectrum of the quantum amplitude: (£ 2 ), for e =1, 5, 7.5 
and 8 Transition to chaos appears at e ps 8. Initial conditions as in Fig. 10. 

Fig. 15: Time evolution of two functions: the unperturbed Hamiltonian of quantum 
linear oscillator (4.3), (Hlo), and the amplitude (£ 2 ), for the values of the perturbation 
parameter e = 0, 0.5 and 2. The corresponding classical dynamics is regular in this 
case, e =0, 0.5, 2. Initial conditions as in Fig. 10. 

Fig. 16: The time evolution of the same functions shown in Fig. 14 but for larger 
values of the parameter e =3, 5, 7.5 and 8. The curves (3) and (4) correspond to the 
classically chaotic regime of motion. Initial conditions as in Fig. 10. 



24 



Berman et al., Fig. 2 



Upper 
\L) ? Manifold 



T 



Detuning 



Laser 
(frequency (O a ) 



M) > 



Lower 
Manifold 



Berman et al Fig. 3 




20.0 



10.0 



0.0 



-10.0 



Phase Space (£, d^/dt); jll=3.99 



8=2.0 



-.1 



20.0 



10.0 



0.0 



-10.0 



8=2.5 



-20.0 l^^^^^^^j _ 20 .o 

-20.0 -10.0 0.0 10.0 20.0 -20.0 -10.0 0.0 10.0 20.0 



8=3.0 



8=4.0 



20.0 



10.0 



0.0 



-10.0 



-20.0 



20.0 



10.0 



0.0 



-10.0 



-20.0 



-20.0 -10.0 0.0 10.0 20.0 -20.0 -10.0 0.0 10.0 20.0 



Phase Space (£, d^/dt); jll=3.99 



8=5.0 



8=8.0 



20.0 



10.0 



0.0 



-10.0 



-20.0 
-20.0 



40.0 
30.0 
20.0 
10.0 
0.0 
-10.0 
-20.0 
-30.0 



-40.0 




-25.0 -10.0 5.0 20.0 35.0 



8=20.0 



-40.0 -25.0 -10.0 5.0 20.0 35.0 



-40.0 -20.0 



0.0 



20.0 40.0 



-2.5 



-10.0 



0.0 



1.0 
0.5 
0.0 
-0.5 
-1.0 



Evolution 5(x); |li=3.99 

8=0.0 



8=0.5 




0.0 20.0 40.0 60.0 80.0 100.0 



8=1.0 




0.0 20.0 40.0 60.0 80.0 100.0 



8=2.0 





-5.0 



-10.0 

20.0 40.0 60.0 80.0 100.0 0.0 20.0 40.0 60.0 80.0 100.0 



8=5.0 



8=8.0 




) i < ^ < ^ < ^ < ^ < 1 _ 1-0 

0.0 20.0 40.0 60.0 80.0 100.0 



20,0 40 60.0 80.0 100.0 0.0 20.0 40.0 60.0 80.0 100.0 

x x 



Fourier Transformation: Re F[^(x)](v); |li=3.99 



8=0.0 



3.0 
2.0 
1.0 
0.0 



-1.0 



0.00 0.50 1.00 

8=1.0 



7.0 



2.0 



-3.0 



0.0 0.5 1.0 

8=5.0 



0.4 
0.2 
0.0 
-0.2 




0.0 0.5 



1.0 

v 



- 7.0 

2.0 

-3.0 

1.50 2.00 0.0 

7.0 



2.0 



-3.0 



1.5 2.0 0.0 



0.1 

0.0 

^^^^ -0.1 
1.5 2.0 0.0 0.5 



8=0.5 



0.5 1.0 

8=2.0 



0.5 1.0 

8=8.0 



1.0 

v 



1.5 2.0 



1.5 2.0 







^—v . A : 







1.5 2.0 



(0):e=0.0; (l):e=0.5; (2):e=1.0; (3):e=1.5; (4):e=2.0 (|u=3.99) 



ic (x)r 




0.40 
0.30 
0.20 
0.10 



A 



(2) 



(4) 



A 



(3) 




\ 



i 1 n-i 



0.030 
0.020 
0.010 
0.000 

( 

0.015 
0.010 
0.005 
0.000 



|c 2 (x)l 



(3) 



(4) 



(2) 



A 



A 1 j 



l\, 



'I 



0.0 



5.0 



10.0 



^ 0.00 
15.0 0.0 



0.020 
0.015 
0.010 
0.005 
0.000 



/ \ A". / i 
A/ '•..-aVA 'K'/,-,/ y 




5.0 



10.0 



15.0 



|c 3 (x)l 




(4) 



\S \\ 



• ^>--A,,^a^' 'A'A^AA 



2^ 



15.0 



0.0 



5.0 



10.0 



15.0 



|c 4 (x)l 



/r* A 




>A/ v ~' >A / - 

A. v"/ 



0.0 



5.0 



10.0 



15.0 



0.010 
0.008 
0.006 
0.004 
0.002 
0.000 



ic 5 (x)r 




0.0 



5.0 



10.0 



15.0 



2.0 



1.5 



1.0 



0.5 



Evolution: 



8=1.0 



o.o 



o.o 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 



0.0 



10.0 



20.0 



8=7.5 



10.0 20.0 

x 



30.0 




30.0 




1.0 



Fourier Transformation: Re F[|c (i)| ](v) 
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Fourier Transformation: Re F[<^ (t)>](v) 
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